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Abstract. Deep photometric data in the V-, R-, I-, z- and K-bands for the cluster 
of galaxies EIS 0048-2942 are used to investigate the properties of the galaxy 
populations at z ~ 0.64 in a field of 2.5 x 2.5 Mpc^. The sample of candidate 
cluster members (N = 171) is selected by the photometric redshift technique and is 
complete up to I = 22.5. Galaxies were classified as spheroids and disks according 
to the shape of the light profile in the I-band, as parametrized by the Sersic index. 

In both optical and NIR, spheroids define a sharp colour-magnitude sequence, 
whose slope and zero points are consistent with a high formation redshift (zf > 2). 

The disk population occupies a different region in the colour-magnitude diagram, 
having bluer colours with respect to the red sequence. Interestingly, we find some 
level of mixing between the properties of the two classes: some disks lie on the 
colour-magnitude sequence or are redder, while some spheroids turn out to be 
bluer. The spatial distribution of cluster galaxies show a clumpy structure, with 
a main over-density of radius ~ 0.5 Mpc, and at least two other clumps distant 
~ 1 Mpc from the center. The various sub-structures are mostly populated by the 
red galaxies, while the blue population has an almost uniform distribution. The 
fraction of blue galaxies in EIS 0048-2942 is fn = 0.11 db 0.07. This is much lower 
than what expected on the basis of the Butcher-Oemler effect at lower redshifts. 

Key words. Galaxies: clusters: individual: EIS 0048-2942 - Galaxies: evolution - 
Galaxies: fundamental parameters - Galaxies: photometry 
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1. Introduction 


Since 1978 (Butcher & Oemler |l978ai |l978bD , several studies have been fruitfully under¬ 
taken in order to investigate the galaxy evolution in clusters. 

Photometric and spectroscopic works found a systematic enhancement of the frac¬ 
tion of blue galaxies in the cores of distant clusters up to z ~ 0.5, known as Butcher 


- Oemler (BO) effect (e.g. Butcher & Oemler 1984, Couch & Sharpies 1987, Rakos & 
Schombert 1995| , Margoniner et al. 2001, Ellingson et al. ^00 1| ), suggesting a strong in¬ 
crease in the star formation of cluster galaxies with redshift (but see also Koo et al. 


It was also found that the blue galaxy fraction increases with cluster-centric distance 


(e.g. Kodama & Bower pOOll) . Recent studies pointed out, however, that selection cri¬ 
teria, as for instance the use of optically selected samples, could affect significantly the 


determination of the BO effect (e.g. Andreon & Ettori 1999 , Eairley et al. 2002 ). For what 
concerns early-type galaxies, photometric studies at different redshifts, showed that they 
appear to be a passively evolving population, where the reddening of more massive galax¬ 
ies with magnitude (the colour-magnitude relation) is the result of a mass-metallicity 
relation (e.g. Aragoh-Salamanca et al. 1993|; Kodama & Arimoto 1997; Gladders et al. 


1998; Stanford et al. 199S; Kodama et al. 1998, hereafter KABA98). 

Different studies have attempted to explain the observational scenario. Kodama & 


Bower (2001) showed that the colour-magnitude (hereafter CM) diagrams of cluster 


galaxies at intermediate redshift can be explained by the continuous accretion of field 
galaxies, whose star formation is inhibited by the hostile cluster environment. The origin 
of the blue galaxies has also been investigated by using spectroscopic and morphological 
data. Barger et al. 


and Poggianti & Barbaro (|1996D showed that a significant frac¬ 
tion of galaxies in distant clusters could have experienced a recent burst of star formation 
(but see also Abraham et al. 1996). Moreover, Couch et al. 


found that most of the 
blue galaxies have the same properties of normal star forming held spirals, while Dressier 


et al. (|1997D showed that the morphological content of distant clusters evolves signih- 
cantly, the fraction of SOs decreasing with redshift. Up to date, it is unknown the origin 
of the mechanism that drives changes of star formation and morphology when galaxies 
are accreted into the cluster, although it is clear the effectiveness of detailed photometric 
and spectroscopic studies in order to explain the scenario of galaxy evolution. 


The cluster EIS 0048-2942 (hereafter EIS0048) was detected by Olsen et al. (1999a) 
using the I-band galaxy catalogue of the ESO Imaging Survey (EIS, Renzini & da Costa 
1997) and the cluster redshift was estimated to be z ~ 0.6 by the matched filter method 


see Olsen et al. |l999b[) . This candidate was also identihed in the EIS catalogue by Lobo 


et al. (|2000|) with an estimated redshift z 0.85 by a modified version of the matched 
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Fig. 1. V-band FORS2 image of the cluster field. The rectangles show the K-band 
ISAAC pointings, while the spatial scale is reported in the lower-left corner. 


filter algorithm. Later, a spectroscopic follow-up (Serote-Roos et al. 2001 ) confirmed the 
presence of a structure at z ~ 0.64 with 15 concordant redshifts. 

This work is the first of a series devoted to the detailed study of this new galaxy 
cluster at z ~ 0.64, for which multi-band photometry and spectroscopic data have been 
collected. 

The paper consists of two parts. The first part deals with data reduction and sample 
selection, while in the second part we discuss the photometric properties of the galaxy 
population. The reader interested only in the second part may skip to Sect. 7. In Sect. 2 
we present the photometric observations. Data reduction and photometric calibration 
are described in Sects. 3 and 4 for VRIz and K-band, respectively. Aperture photometry 
is presented in Sect. 5, while Sect. 6 deals with the photometric redshift estimates. In 
Sect. 7 we introduce the galaxy shape classification. The colour-magnitude distributions 
are discussed in Sect. 8. The spatial distribution of the cluster population and the BO 
effect are studied in Sects. 9 and 10, respectively. A summary is given in Sect. 11. 

In the following we assume ftm = 0.3, Da = 0.7 and Hq = 70 Km s“^Mpc“^. With 
this cosmology the age of the universe is ^ 13.5 Gyr, and the redshift z = 0.64 corre¬ 
sponds to a look-back time of ~ 6 Gyr. 


2. Observations 

New photometric observations of the cluster of galaxies EIS 0048 were carried out at the 
ESO Very Large Telescope (VLT) during two observing runs in August 2001. All the 
nights were photometric with excellent seeing conditions. 

The data include VRIz and K-band imaging taken with FORS2 and ISAAG instru¬ 
ments, respectively. The VRIz images consist of a single pointing of 6.8 x 6.8 arcmin^ 
for each band, while for the K-band a mosaic of four pointings covers a total area of 
4.9x4.9 arcmin^. Furthermore, we obtained a mosaic of three FORS2 pointings in I-band 
by using the high resolution (HR) observing mode. The reduction of the HR images will 
be presented elsewhere. In order to correct for the fringing patterns in the I(SR)- and 
z-bands, we obtained five dithered exposures with a dithering box of size 25". For the 
K-band, we collected two sets of 22 exposures for each pointing with DIT = 12 sec and 
NDIT = 6, by using a dithering box of size 15". The overlapping of the VRIz and K-band 
images is shown in Fig. while the relevant information on the data are summarized in 
Tabled 
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Table 1. Relevant information on the observations. Col. 1: waveband. The symbols SR 
and HR denote the standard and high resolution I-band images, respectively. Col. 2: pixel 
scale. Col. 3: total exposure time for each pointing. Col. 4: full width half maximum of the 
seeing disk. For the I(HR)- and K-bands, we report the range for the various pointings. 


Band 

Sc9rl6 Texp 

"/pyi\ ksec 

Seeing 

n 

V 

0.2 

3.6 

0.6 

R 

0.2 

2.4 

0.9 

I(SR) 

0.2 

1.2 

0.8 

I(HR) 

0.1 

1.8 

0.3-0.4 

z 

0.2 

1.8 

0.6 

K 

0.15 

3.2 

0.4-0.5 


Fig. 2. Correction of the fringing in the z-band image. Upper - left and right panels show 
a portion of the z-band image before and after the correction. The profiles of the average 
of three columns are shown in the lower plots. Intensities are in units of the standard 
deviation of the background, ctbg- 

3. VRIz data 

3.1. Reduction 

The data reduction was performed by using IRAF ^ and Fortran routines developed by 
the authors. 

The bias was corrected by using zero exposure images and the prescan/overscan 
regions of the frames. For each filter, the scientific images were divided by a flat-held 
frame obtained by combining twilight sky exposures. After this correction, however, low 
frequency variations of ~ 5% were still present across the images. We obtained, therefore, 
a sky flat frame by median combining the scientific exposures taken during the night. 
In these frames, residual objects were masked and a polynomial fit was applied. The 
division of the scientific frames by the fitting surfaces allowed to reduce the low frequency 
variations across the chip to ^ 0.5%. 

After this procedure, fringing patterns were noticeable in the R-, T and particularly 
in the z-band frames. For the R-band, we did not apply any correction since the fringes, 
found only in the right side of the chip, showed peak-to-peak variations in the final image 
^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated 
by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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of less than half the r.m.s. of the background. For the I-band, fringing was almost com¬ 
pletely removed because of the dithering between the various exposures. Since dithering 
did not remove the fringing pattern in the z-band, we created a fringing model by hltering 
the relative sky flat frame. The model was suitably scaled and then was subtracted from 
each scientific image, as illustrated in Fig. |^. This allowed to reduce the peak to peak 
variations of the fringing pattern in the final image from ~ 4 times to ^ 0.5 times the 
background noise. 

For each image, cosmic rays were detected by using the IRAF task COSMICRAYS 
and a mask frame was produced including cosmic rays and chip defects. The images 
were corrected for the different airmass with the relative extinction coefficients (see next 
section) and added by using the corresponding masks with the IRAF task IMCOMBINE. 


3.2. Photometric calibration 


The VRI images were calibrated into the Johnson-Kron-Cousins photometric system 


by using comparison standard fields from Landolt ( 1992 ) observed during each night. 
Instrumental magnitudes of the standard stars were derived within an aperture of 10'' 


by running SExtractor (Berlin & Arnouts |1996|) . Eor each waveband, we adopted the 
following calibration relation: 


M = M'-b 7 -C-A-X-tZP 


( 1 ) 


where M and C are magnitude and colour of the standard star, M' is the instrumental 
magnitude, 7 is the coefficient of the colour term, A is the extinction coefficient, X is 
the airmass, ZP is the zero point. For each waveband, the quantities 7 , A and ZP were 
derived by a robust least square fit to the data, except when standard stars with a 
sufficient range of colours and/or airmasses were not available. In such case, we adopted 
the value of 7 and/or A of the FORS2 standard calibrations!^ relative to the period of our 
observations. In Table we report the fitting coefficients, the corresponding values of the 
FORS2 standard calibrations and the uncertainty on the photometric calibration given 
by the r.m.s. of the residuals to the fit. The results of the fit are in very good agreement 
with the FORS2 standard calibrations, with differences in the zero points smaller than 
0.03 mag. 

To check the accuracy of the photometric calibration, we considered the distribution in 
the (V-R,R-I) diagram of the stars present in the cluster field. The observed colours were 
corrected for the different seeing of the images as described in Sect. 5. In Fig. ^ (upper 
panel) we compare this distribution with that obtained from the synthetic stellar spectra 
of the HILIB library (Pickles |1998D . The synthetic colours were produced by using the 
total efficiency curves of FORS20 and were calibrated with respect to the Vega spectrum. 


^ http: / / www.eso.org/observing/dfo/quality/FORS/qc/zeropoints/zeropoints.html 
® See the VLT exposure time calculator at http://www.eso.org/observing/etc/ 
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Fig. 3. Colour - colour diagrams of the stars in the cluster field (black circles). Gray 
circles mark the synthetic colours for the stellar spectra of Pickles (|l998|) . The curve in 
the lower panel is the polynomial fit to the locus of the Pickles stars, that was used to 
match the observed and the synthetic distributions. 


The observed and the synthetic distribution are in very good agreement, confirming the 
accuracy of data reduction and calibration, and the photometric quality of the observing 
nights. 

Since no standard comparison fields were available for the z-band, we adopted the 
following calibration method. Synthetic z magnitudes were produced from the Pickles 
spectra and were calibrated by assuming I — z = 0 for the Vega spectrum. The zero point 
was obtained by fitting the distribution in the (R-I,Tz) diagram of the stars in the cluster 
field, corrected for seeing effects, to the locus of the synthetic colours. The result of the 
fit is illustrated in Fig. (lower panel). The uncertainty on the zero point was estimated 
with the bootstrap method by taking into account the uncertainties on the colours of 
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Table 2. Results of the photometric calibration of the VRI data. Col. 1: waveband. Col. 2: 
colonr C used in Eq. (^). Cols. 3, 4: colour term coefficients 7 and 7 '. Cols. 5, 6 : atmospheric 
extinction coefficients A and A'. Cols. 7, 8 : zero points ZP and ZP' scaled to Isec. Col. 9: r.m.s. 
of the fitting residuals. The quantities marked by a prime are the FORS2 standard calibrations 
(see note ^). The uncertainties refer to Icr standard intervals. Dots denote the quantities kept 
fixed in the fits (see text). 


C 7 7 ' A A' ZP ZP' (7 

V V-I . 0.010 ±0.005 0.18 ±0.02 0.154 ±0.007 27.280 ± 0.030 27.287 ±0.014 0.005 

R V-R . -0.018 ±0.011 0.10 ±0.01 0.114 ±0.007 27.420 ± 0.013 27.450 ± 0.015 0.013 


I V-I -0.058 ±0.004 -0.059 ±0.005 . 0.053 ± 0.007 26.474 ± 0.003 26.484 ± 0.037 0.018 


the stars. We obtained ZP^ = 25.216 ± 0.012. For the atmospheric extinction coefficient 
in the z-band, we assumed A^ = 0.08 mag/airmass, that is typical for the Paranal site. 


4. K-band data 

4.1. Reduction 

The data reduction was performed by using ISAAC pipeline recipes included in Eclipse, 


version 4.0, the InfraRed Data Reduction (IRDR) software (see Sabbey et al. 2001 ) and 
Fortran routines developed by the authors. 

For both observing nights, we derived a dark frame by using the dark recipe of Eclipse. 
Since the dark variations between the two nights were less than 1 — 2%, the dark frames 
were averaged and then subtracted from all images. As shown in Fig. a double ramp 
pattern along the column direction was noticeable in the scientific images after dark 
subtraction. The effect was found to be almost constant along the chip rows, varying 
from ~ —5% at the rows where the detector starts to be read out (rows 1 and 513), to 
~ 3 — 5% at the rows read last. Since this pattern was not present in the flat-field frame, 
it must be a residual additive signal uncorrected by the dark subtraction. As discussed 
in the ISAAC Data Reduction Manual (version 1.5, Amico et al. 2002| ), this effect can be 
due to the dependece of the ISAAC infrared detector bias on time and on illumination 
flux. To remove such an effect, we applied a suitable correction immediately after dark 
subtraction. Each image was collapsed by taking the median along the row direction, and 
a low-pass filter was applied to the first 512 and to the last 512 samples of this signal 
in order to remove high frequency variations due to the noise and to the pixel to pixel 
response of the chip. The one-dimensional frame obtained by this procedure (see Fig. ^ 
was subtracted column by column from the original image, removing the ramp pattern 
at better than ~ 1%. 
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Fig. 4. Pattern remaining in the K-band images after dark subtraction. The gray line 
represents the signal after low-pass filtering (see text). 


For the flat-field correction, we divided the scientific images by a differential flat- 
field frame, obtained by two sets of twilight sky exposures taken in the morning and in 
the evening, respectively. Each twilight exposure was corrected for the odd-even effect, 
particularly noticeable in the lower left quadrant (where it amounted to ~ 2% peak to 
peak variation). To this aim, we used the oddeven recipe of Eclipse. For each set of twilight 
images, a flat-field frame was obtained by subtracting high and low counts exposures in 
order to correct for residual additive components not removed by the dark subtraction. 
The final flat-held was derived by averaging the morning and evening twilight flats. This 
procedure allows to minimize intensity variation within the array at twilight. To test the 
accuracy of the flat-field correction, we compared the magnitudes of each photometric 
standard star at the different chip positions. The differences were found to be ~ 5 — 7% 
(Icr standard level). To achieve higher accuracy, we applied an illumination correction 
by retrieving from the ESO archive the illumination frames closest in time to the date 
of our observations. After illumination correction, the variation were found to be ^ 1%, 
that we assumed as the final accuracy of the flat-field. 

The subsequent reduction included sky subtraction and image combining, that were 
performed by a two-step procedure. First, the images were processed by using the IRDR 
software. For each frame, a sky image was obtained by a robust mean of the eight nearest 
sequence exposures. After sky subtraction, a mask image was created for each frame by 
running SExtractor (Bertin & Arnouts 1996) with the checkimage OBJECTS option 
and was used to measure the dithering offsets. The first step coadded images were then 
obtained by combining the exposures within each sequence. At the second step, we created 
an object mask for the coadded images obtained at the first iteration by using SExtractor. 
These masks were expanded by a factor of 1.5 and de-registered to the corresponding 
dither exposures. Pixels with low and high counts in the flat-field were also included 
into the masks in order to reject hot and cold pixels. For each frame, the sky image was 
obtained by the average of the six nearest exposures rejecting masked pixels. The images 
were then coadded by using the IRAF task IMCOMBINE with a SIGCLIP algorithm for 
cosmic ray rejection. 


4.2. Photometric calibration 


The photometric calibration of the ISAAC data was performed by using standard stars 
from the list of Infrared NICMOS Standard Stars (Persson et al. 


and from the list 


of UKIRT standard stars (see Hunt et al. 1998). 
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A(l) 


Fig. 5. Comparison of I-band high resolution and standard resolution aperture magni¬ 
tudes. Differences are in the sense SR - HR. The dot - dashed line marks the aperture 
correction estimated by the method described in the text. 


The images were airmass-corrected by assuming an extinction coefficient 
Ak = 0.08 mag/airmass, typical for the period of our observations, and the total fluxes 
of the standard stars were estimated within a circular aperture of diameter 8". A small 
systematic difference of ~ 0.05 mag was found between the zero points derived from the 
UKIRT and NICMOS standard stars. Such result is not unexpected, since the UKIRT 
standards can be affected by small systematic and/or random errors with respect to 


those of the NICMOS list (see Sect. 5.2 of Persson et al. 1998 ). The final zero point was 
obtained, therefore, by considering only the NICMOS standard stars, and was estimated 
to be 24.34 ± 0.01 mag (scaled to Isec exposure time). This is in good agreement with 
the value of the ISAAC pipeline (ZP ~ 24.36), relative to the period of our observations. 


5. Aperture photometry 

A catalogue of each image was produced by using the software SExtractor (Bertin & 


Arnouts 1996). Since the images have different seeing, particular care was taken in the 
choice of deblending parameters in order that objects were deblended in the same way for 
each band. For each object, we obtained magnitudes within a fixed aperture of diameter 
2.3", corresponding to a physical size of ^ 15 kpc at z = 0.64, and magnitudes within an 


aperture of diameter a ■ rx, where rx is the Kron radius (Kron 198C). We chose a = 2.0, 


for which the Kron magnitude is expected to enclose 92% of the total flux. The total 
magnitudes were computed by adding 0.08 mag to the Kron magnitudes. 
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The catalogues were cross-correlated by deriving the coordinate transformation be¬ 
tween each image and the V-band image. To this aim, we proceeded as follows. A first 
list of matched objects was obtained by using the IRAF task XYXYMATCH, with the 
triangle pattern algorithm. This method has the advantage to be robust with respect to 
geometric distortion effects in the images. The final transformations were obtained by 
fitting the coordinates of the objects in the matching list with polynomials of order 5, 
which allowed to correct for geometric distortion effects with an accuracy better than 
0.5 pixel. The derived coordinate transformations were used to match the catalogues by 
using a nearest algorithm, producing a final list of N = 4868 sources detected in at least 
one band. For the objects with multiple K-band observations, the magnitudes turned 
out to be fully consistent within the photometric errors, and the final magnitudes were 
estimated by a weighted mean of the various measurements. 

Since the images have different seeing (cfr. Table |^), the derivation of the colour 
indices requires suitable aperture corrections. To this aim, each image was convolved 
with a double gaussian kernel, whose parameters were chosen interactively to match 
the PSF of the image with the worst seeing (R band). Noise was added in order to 
compensate for the smoothing introduced by the convolution. The aperture corrections 
were estimated by the robust mean of the differences between the magnitudes of the 
objects in the original and the convolved images, and amount to 0.05 mag in V, 0.03 mag 
in I, 0.06 mag in z, and 0.10 mag in K. The corresponding uncertainties were estimated 
by the errors on the mean, and are 0.007 mag, 0.013 mag, 0.02 mag, and 0.03 mag in V, 
1, z and K, respectively. These values were added in quadrature to the photometric errors 
of the magnitudes. To test our procedure, we compared the magnitudes of the sources in 
common between the standard and high resolution 1-band images, whose seeing values 
are 0.3" and 0.8", respectively. The magnitude differences are shown in Fig. together 
with the aperture correction derived as described above. The figure demonstrates the 
accuracy of the adopted method. 


The completeness in each band was estimated following the method of Garilli et 
al. ( 1999| ), that consists in the determination of the magnitude at which the objects start 
to be lost since they are below the brightness threshold in the detection cell. Details will 
be given elsewhere (Massarotti et al. 2002, in preparation). The completeness magnitudes 
in V-, R-, I-, z- and K-band are 25.7, 25.0, 23.2, 22.5 and 21.2 respectively. 


6. Photometric redshifts 

Photometric redshifts were estimated according to the Spectral Energy Distribution 
(SED) fitting method (see Massarotti et al. 2001a, b, and references therein). In or¬ 
der to achieve a reasonable accuracy, we considered galaxies with signal-to-noise ratio 
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Fig. 6. Distribution of photometric redshifts Zp for galaxies with 18.5 < I < 22.5, after 
subtraction of the field counts estimated by the VIRMOS survey data. The dashed line 
marks the redshift 0.64. 


S/N > 5 in at least three bands. This sample contains N = 633 galaxies, of which N = 360 
have K-band information, and is complete up to the total magnitude It = 22.5. 

The photometric redshifts were derived as described in Busarello et al. ( ^002 ), here we 
briefly outline the technique. Taking into account the depth of our imaging, we looked for 
redshifts in the range z G [0.0, 2.5] with a step of 0.01. Model galaxy spectra were provided 
by the code of Bruzual & Chariot ( |1993 ). The adopted templates consist of models with 
a Scalo (1986) IMF and with an exponential SFR . We chose t = 1,4,15 Gyr, to 
describe the colours of EjSQ, Sa/Sb, and Sc/Sd templates. To allow for different metal- 
licities of early-type galaxies, we introduced E/SO models with Z/Zq =0.2, 0.4, 1 and 
2.5, while template spectra evolution was followed in the time interval t G [0.01,12.0] 
Gyr. To estimate the uncertainty Azp in the photometric redshift, we performed nu¬ 
merical simulations by taking into account the uncertainties on magnitudes. Since the 
K-band does not cover the whole cluster field, we investigated the possible bias on Zp for 
the galaxies without K-band data. To this aim, we compared the photometric redshifts 
obtained for the sources with It < 22.5 (a) by taking into account and (b) by ignoring 
the K-band. The number of galaxies which were found to be cluster members 
in case (a) and no cluster members in case (b), and viceversa, turned out to 
be negligible (< 5%). 

We defined a galaxy as a cluster member when the photometric redshift is in the 
range Zp G [0.54,0.74]. However, to obtain a reliable sample of cluster members via pho¬ 
tometric redshifts it is crucial to take into account possible degeneracies of template 
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colours at different redshifts. By inspecting colours of GISSEL98 templates, we found 
that some degeneracy exists between colours of blue Sc/Sd galaxies at the cluster redshift 
(V — R = 0.94, R — I = 0.94,1 — K = 2.62 with tform = 12 Gyr, Z = Zq) and those of red 
if/5'0 galaxies at z ~ 0.3 (V — R = 1.04, R — I = 0.82,1 — K = 2.80 with tform = 12 Gyr, 
Z = Zq). To correct for this effect, we also considered as cluster members the objects 
having z G [0.3,0.4] and Azp consistent with the cluster redshift. We found that N = 19 
galaxies satisfy this criterion. Most of these galaxies (~ 75%) turn out to be disk domi¬ 
nated objects and therefore their late spectral type is also supported by the shape of the 
light profile (Sect. 7). The final list of cluster members brighter than I = 22.5 consists in 
N = 171 galaxies. 


To gain insight into the contamination by field galaxies in our sample, we used the 
VIRMOS preparatory photometric survey as the control sample. These data, kindly pro¬ 
vided by the VIRMOS Gonsortium (Le Fevre et al. 2002, in preparation), are actually 
the proper dataset for estimating the field contamination in our sample because of their 
large field, similar filters and high depth). In particular, to estimate the redshift distribu¬ 
tion for field galaxies we used the VIRMOS catalogues (VRIK, McGracken et al. 2002, in 
preparation) relative to a 150 arcmin^ field with available K-band photometry. It has to 
be noticed that the photometric redshifts were obtained for the field and the cluster by 
using a similar photometric baseline and the same procedure. The distribution of photo¬ 
metric redshifts, after subtraction of the field counts, is shown in Fig. for the objects 
with 18.5 < I < 22.5, corresponding to the magnitude range of cluster galaxies up to the 
completeness limit of the sample with photometric redshifts. It is dominated by the peak 
around the cluster redshift. The FWHM of the peak is ^ 0.1 and is due to the intrinsic 
redshift distribution of the cluster members and to the errors on photometric redshifts. 
The redshift distribution of VIRMOS galaxies is in full agreement with that relative to 
the cluster field outside of the peak at z ^ 0.64, thus confirming the effectiveness of 
our selection procedure and the reliability of the VIRMOS control field. In the field of 
FIS 0048 we found N = 140 ± 12 galaxies with z < 0.5 and N = 46 ± 7 galaxies with 
0.8 < z < 1.8, while in the same redshift ranges the number of galaxies predicted on the 
basis of the control sample are N = 141.5 ± 6.6 and N = 43 ± 4 respectively. According 
to the VIRMOS data, N = 48 ± 4 (25%) galaxies of the cluster sample are expected to 
be field galaxies. 


In Fig. ^ we compare photometric (zp) and spectroscopic (zs) redshifts for a sub¬ 
sample of 25 cluster galaxies (Busarello et al. 2002, in preparation, and C. Lobo, private 
communication). We notice that all the spectroscopically confirmed cluster members 
satisfy the selection criterion previously discussed. 
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Fig. 7. Comparison of spectroscopic Zg and photometric redshifts Zp for 25 spectroscop¬ 
ically confirmed cluster members. 


7. Galaxy shape classification 

Although a proper morphological classification requires the high quality of the HST 
images, relevant information on the properties of the light distribution of distant galaxies 


can be obtained by using ground based data (see La Barbera et al. 2002 , hereafter LBM02, 
and references therein). We used the procedure described in LBM02 to model the surface 
brightness distribution of galaxies in the cluster field, by fitting Sersic models convolved 
with the PSF of our images. This analysis will be presented in detail in a forthcoming 
paper. 

For the present work, we will use the information on the shape of the galaxy light 
profiles as parametrized by the Sersic index n. We fitted Sersic models both to the stan¬ 
dard resolution (SR) and to the high resolution (HR) I-band images. Accurate structural 
parameters for the galaxies at the redshift of EIS0048 can be obtained from the very 
good seeing of the HR images (La Barbera et al. 2002, in preparation). Nevertheless, we 
found that a reliable separation between galaxies with low n value (disks) and galaxies 
with a higher value of the shape parameter (spheroids) can be also done with our SR 
images, allowing us to perform the galaxy shape classification in the whole observed field. 
In fact, only ^ 8% of the galaxies brighter than It = 22 have a different classification 
in the HR and the SR images, while the number of mismatches increases to ~ 17% at 
It ^ 22.5. Moreover, for the galaxies with discordant classification we did not find any 
significant systematic effect: ^ 10% of the galaxies defined as disks on the basis of the 
HR data are classified as spheroids in SR, and vice-versa for the remaining ~ 7%. 
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Fig. 8. Colour-magnitude diagrams for the galaxies in the cluster field (gray circles). 
Black symbols denote objects with It < 22.5 defined as cluster members on the basis 
of their photometric redshift. Circles denote the spheroids, while disks are marked by 
crosses. The dashed lines are the best fits to the CM while the dotted lines limit the ±cr 
intervals. The error bars in the left corners of the plots indicate the mean uncertainties 
on colours and magnitudes at It ~ 19.7 and It ~ 22.2 for the upper plot and at Kt ~ 17 
and Kt 19 for the lower plot. 


We chose n = 2 to discriminate between the two classes of galaxies. As shown by 


numerical simulations (see van Dokkum et al. 1998), this criterion corresponds to separate 
objects with a low bulge fraction (< 20%) from galaxies with a larger bulge component. 
Our final sample of galaxies with shape information consists of all the N = 171 cluster 
members brighter than It = 22.5. 


8. Colour-magnitude relations 

The optical (V-I) and NIR (V-K) colour-magnitude diagrams for the objects in the cluster 
field are shown in Fig. |[ The black dots denote objects brighter than It = 22.5 defined 
as cluster members according to their photometric redshift (see Sect. ||), while different 
symbols are used to mark disks and spheroids. At this magnitude limit, we have shape 
classification for all the cluster members. 


8.1. Distribution in the CM planes 

Spheroids define sharp CM relations both in the (V-K,K) and in the (V-I,I) planes. 
The sequences extend for at least 3 magnitudes, up to the completeness limits of our 
samples, and are defined in the plots by the loci enclosed within the dotted lines (see 
next section for details). It is interesting to notice that 43 spheroids are significantly 
bluer with respect to the V — I red sequence, among them 24 are bluer in the range 
[—0.4, —0.2] mag. In particular, we find six blue spheroids at very bright magnitudes 
(It < 20.2 and 1.9 < V — I < 2.2). Two of these have also K-band photometry and lie 
on the red locus in the (V-K,K) plane. We verified that the classification of these blue 
spheroids is the same in all our images, and therefore it is not affected by a possible 
uncertainty on the Sersic index. Most of the disks (43) are bluer with respect to the 
optical red locus, although 10 galaxies do not show significant deviations and other two 
have redder colours. Four of these 12 galaxies have K-band photometry and are located 
above the V-K CM relation. 

Since our sample is selected by the photometric redshift technique, some consider¬ 
ations are needed about the contamination by field objects in our diagrams. On the 
VIRMOS data, in the (V-I,I) plane we expect 39.5 ± 3.5 blue galaxies, 8 ± 1.6 red se- 
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Fig. 9. I-band image of the cluster field. The map of the counts density is represented 
by the gray contours. The grey intensity is proportional to the density value. Triangles 
and circles mark disks and spheroids respectively. 

quence objects and 0.6 ±0.4 redder galaxies. These numbers indicate that the presence of 
blue spheroids and red disks is not due to the field contamination but that at least some 
of these objects are real cluster members. Moreover, on the basis of the comparison field 
sample, the cluster membership of the bright blue spheroids is highly significant. In fact, 
in the range of magnitude and colours of these galaxies, we expect that only 0.5 ± 0.5 
objects could be field contaminants. 


8.2. Slope and zeropoint of the CM relations 

We study the CM relation of EIS0048 by comparing our results with the models of 
KABA98. We performed a linear regression of the V-K and V-I colour-magnitude se¬ 
quences by considering the following relations: 

V-I = Avi ± Bvi • I (2) 

V — K = Avk ± Bvk • K. (3) 

The fits were done by minimizing the bi-weight scatter of the residuals (see Beers, Flynn 


& Gebhardt 1990), and the algorithm was modified in order to weight each point with 
the corresponding uncertainty on the photometric redshift estimate. Only objects within 
the completeness limit of our photometry were considered. This procedure allows to 
minimize the effect of outliers and of objects with poor photometric information and/or 
cluster membership. We also verified that the fitting results do not change by excluding 
disks from our sample. The uncertainties on the fitting parameters were obtained by 
numerical simulations of the distributions in the CM planes. For the optical relation, we 
obtained Ayi = 3.95 ± 0.30 mag and Byi = —0.073 ± 0.015, with a r.m.s. of the residuals 
that amounts to tryi = 0.10 ±0.03 mag. For the NIR, due to the smaller sample size, 
the relative uncertainty on the CM slope is large: Byx = —0.15 ±0.08. Therefore, we 
chose to repeat the fit by using a fixed value of the slope Byx = —0.11, that is typical 
for the spheroids in clusters at the same redshift (see Fig. 4 in KABA98). This gives 
Avk = 7.75 ± 0.05 mag, with r.m.s. of the residuals cvk = 0.14 ± 0.06 mag. The fits to 
the CM relation are shown in Fig. || along with the relative 3 a intervals. We denote an 
object as a blue galaxy if it lies more than 3 a below the (V-1,1) CM fit, as a red galaxy 
otherwise. 

By considering the previous results, we see that the dispersions around the optical 
and NIR CM relations can be fully explained by the typical uncertainty on galaxy colours 
(~ 0.1 mag, cfr. Fig. 0 ), and that, therefore, the intrinsic scatter of the relations is very 
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small. Unfortunately, due to the field contamination and to the uncertainty on morpho¬ 
logical selection, the quoted errors on the dispersions are too large to obtain an accurate 
quantitative estimate of the intrinsic scatter. The fitting coefficients of the Eqs. (|^) can 
be compared with the predictions of the KABA98 models relative to the V 555 — I 814 vs. 

and the V — K vs CM relations considered by the authors (see Figs. 4, 5 in 
KABA98), where V 555 and I 8 i 4 are the magnitudes of the HST filters. By comparing the 
value of Bvi with the values shown in Fig. 4 (right middle panel) in KABA98, we see 
that our estimate is fully consistent with the prediction of a pure metallicity sequence 
with an old epoch (zf = 4.5) of galaxy formation. Pure metallicity models, with a more 
recent formation epoch, Zf = 1.7 and Zf = 1.2, are discarded at 1.5cr and 2 .Oct confidence 
level respectively. The zero points of the CM relations can be compared with the values 
shown in Fig. 5 in KABA98. To this aim, we referred the Ayi and Avk coefficients to 
the magnitudes Kt = —25.5 and It = —22, by using the same cosmology adopted by 
KABA98 (Ho = 50 Km s“^Mpc“^, = 0.3, Ha = 0.7) and by estimating K correc¬ 


tions from Poggianti (1997). For the NIR, we obtain a corrected value of the zero point 
Ayjj = 5.58 ± 0.05 mag, that is consistent with the models of KABA98 having Zf > 1.7, 
while younger models are not consistent with our data. For the optical zero point, we ob¬ 
tain Ayj = 2.33 ± 0.03 mag that is about 0.2 mag below the value relative to the model 
with Zf ~ 4.5. We notice, however, that this difference can be affected by the uncertainty 


on the zero point of the HST photometry (see e.g. Holtzman et al. 1995), making this 
disagreement not significant. In fact, the comparison of the R — K CM with KABA98 
gives the same results obtained for the V — K relation. Moreover, we notice that the 
V — I colour predicted by the GISSEL98 synthesis code (see Bruzual & Chariot 1993| ) for 
an E template, with solar metallicity and with Zf = 4.5, is 2.34 mag, in full agreement 


with the value of A- 


VI- 


9. Spatial distribution of cluster populations 

To further investigate the properties of the blue and red galaxy populations, we analyzed 
their spatial distribution in the cluster field. In Fig. |^, we show a map of the number 
density of the objects defined as cluster members according to their photometric redshifts. 
Cluster galaxies show an evident clumpy distribution: the peak of the density is associated 
to a main structure having projected radius ^ 1.0' — 1.5' (^ 0.5 Mpc at z = 0.64), while 
at least two secondary clumps of galaxies are found at a distance of ~ 3.0' (~ 1.2 Mpc) 
from the main peak. Although we cannot quantify the number of spheroid and disks that 
are expected to be field objects, it is quite remarkable that the various structures are 
populated, mostly, by galaxies with a more concentrated shape, while disks are found 
preferably in the low density regions. 
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Fig. 10. Map of the counts density for the red and blue cluster populations (upper and 
lower panels). The grey scale is the same for both panels. 

Fig. 11. Density profile of the galaxies defined as cluster members. Error bars denote 
Poissonian uncertainties. Density have been arbitrarily scaled in order to have a density 
of 1 gal/arcmin^ for the field (dashed horizontal line). The curves show the models fitted 
to the density profile (see upper-right). 


By using the separation between red and blue galaxies discussed in Sect. ^ 


we con¬ 


structed the density maps relative to the two classes of objects (see Fig. The clumps 
of galaxies defining the cluster structure are evident in the plot relative to the red pop¬ 
ulation, for which we expect a very low field contamination (~ 7%). With the exception 
of a possible clump near the border of the field, the density map of the blue population 
does not show any particular structure. This is partly due to the field contamination, but 
it could also indicate that cluster blue galaxies do not share the same spatial distribution 
of the red ones. In fact, if we construct the density map of the blue galaxies, by weighing 
each object with the relative luminosity, we find that the clumps present in the lower 
panel of Fig. outside the region corresponding to the cluster central structure, become 
the dominant features in the plot. We point out, however, that a spectroscopic survey in 
the field is needed in order to definitely exclude the contamination effect. 


10. Butcher-Oemler effect 

In order to derive the fraction of blue galaxies fe in the cluster EIS 0048, we followed the 
procedure of Butcher & Oemler ( 1978a , 1978b| ). We computed the fraction of galaxies 
that are 0.2 mag bluer than the red sequence in the rest-frame (B-V,V) diagram within 
a magnitude limit My = —20, and that are enclosed within R 30 , the radius containing 
30% of the total galaxy population. 

In Fig. we show the density profile of the cluster members derived by considering 
circular coronae centered on the peak of the density map (see Fig. ^) with the suitable 
correction for the rectangular geometry of our field. The density counts relative to the 
cluster core decrease gradually up to a distance of ~ 2 arcmin, while the substructures 
of Fig. ^ are shown by the increase of the profile between R = 2' and R = 4'. At large 
distances, the density is very close to the value expected in the field, showing that the 
estimate of field contamination is reliable. To derive the value of R 30 , we need to predict 
the fraction of cluster galaxies fo expected out of the observing field. To this aim, we used 
an extrapolation of the density profile, by fitting the density counts of the cluster core 


(R < 2') by different models. We considered the models adopted by Adami et al. (1998) 
to analyze a large sample of nearby clusters: the generalized Hubble and King profiles, the 


generalized NFW (see Navarro, Frenk & White 1995, 1996) and de Vaucouleurs models. 
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Table 3. Standard parameters for the computation of the BO effect. See text for details. 



c 

AC 

R 30 

AR 30 




(') 

(') 

1(R < 2') 

0.50 

0.02 

1.0 

0.07 

to 

A 

0.43 

0.02 

1.5 

0.17 


The fitted profiles were integrated from R = 4' up to an Abell radius (3h^Q^Mpc), giving 
fo = 20 ± 5%. This value was used to calculate R 30 and the concentration parameter 
C = log(R 6 o/R 2 o)) where Rp is the radius that encloses p% of the total number of cluster 
galaxies. The same procedure was repeated by extrapolating the density profile out of 
R = 2', without considering the substructures present between R = 2' and R = 4'. 
In Table || we report the values of C and R 30 estimated for case 1 (R< 2') and for 
case 2 t(R< 4'), and the corresponding uncertainties, obtained by taking into account 
the Poissonian error on our number counts and the uncertainty on the density of field 
counts. As it was expected, the value of R 30 is smaller for case 1, and ranges between 
R 30 ~ 1' and R 30 ~ 1.5', while the value of C varies from ~ 0.43, that is typical of a not 
particurarly compact cluster, to ^ 0.5, that is closer to that of a more rich, concentrated 


structure (see Butcher & Oemler |l978b[ ). To check our estimate of R 30 and C, we repeated 
the calculation by considering only the red cluster members. We obtained R 30 = 1.2±0.1 
and C = 0.53 ± 0.02 for case 1, while R 30 = 1.5 ± 0.19 and C = 0.44 ± 0.013 for case 2. 
These values are fully consistent with those reported in Table ||. 

The blue fraction fe was computed by using the V — I colour index, that is very 
close to B — V rest-frame at the cluster redshift. In order to convert V — I colour into 
B — V rest-frame, we used the synthesis code of Bruzual & Chariot (1993) to construct 
galaxy templates with different age, metallicity, star formation rate, and dust content, 
and we computed for each of them V — I and B^ — Vr colours, where B^ and Vr are 
magnitudes in the redshifted B and V filters. The relation between the two colours was 
found to be described very well by a polynomial of degree four, with a very small r.m.s. of 
0.025 mag. The same procedure was adopted to convert the I magnitudes into rest-frame 
V-band magnitudes. The value of fe was computed by subtracting the fraction of blue 
galaxies obtained for the cluster sample to that expected for the field. The calculation 
was performed by considering only galaxies brighter than My = — 20 , that corresponds 
to I ^ 22.5 at z = 0.64. We notice that this corresponds to the completeness limit of our 
cluster sample. The uncertainty on fe was estimated by taking into account Poissonian 
and photometric errors of both cluster and field samples. In Fig. |^, we show the trend 
of fe versus the cluster-centric distance. At distances between R = 1' and R = 1.5', that 
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Fig. 12. Fraction of blue galaxies within R, where R is the cluster-centric distance. 


enclose our estimate of R 30 (cfr. Table y), we obtain fe = 0.11 ± 0.07. We notice that this 
result is very robust with respect to our computation of R 30 : fe is, in fact, practically 
constant between R = 1' and R =2'. 

The blue galaxy fraction of the cluster EIS0048 can be compared with the value 
predicted by the standard BO effect at z = 0.5, fe = 0.25, and the values obtained by 
recent works, who have attempted to compute the BO effect at z > 0.5. In particular, 
the BO effect has been investigated at z ~ 0.8 by van Dokkum et al. (2000D, who found 


fe = 0.22 ± 0.05, and by Rakos & Schombert ( 1995 ), who estimated a very high blue 
galaxy fraction, fe ~ 0.8. In agreement with van Dokkum et al. ( ^000 ), our results do 
not indicate a high increase in the BO effect at z > 0.5. The cluster EIS 0048 at z = 0.64 
shows a quite low blue galaxy fraction, about one half of the standard BO effect at z = 0.5. 
By looking at Fig. |^, it is also interesting to notice that the value of fe seems to increase 
at distances greater than R = 2' up to ~ 23% at R = 4'. This would be consistent with 


what found by previous studies at lower redshift (e.g. Kodama & Bower 2001, Fairley et 
al. p 002 ), and can be explained as a consequence of accretion of field galaxies into the 
cluster structure. 


11. Summary. 

We have studied the properties of the galaxy populations in the cluster EIS 0048 at 
z = 0.64, by using a large photometric baseline including V-, R- T, z- and K-band 
data. Cluster members have been selected by using the photometric redshift technique, 
producing a final list of N = 171 galaxies complete up to I = 22.5. To estimate the field 
contamination, we used as control sample the VIRMOS preparatory photometric survey, 
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for which photometric redshifts were estimated by using the same procedure adopted for 
the cluster sample. Out of the N = 171 candidate members, N=48 objects are expected 
to be field contaminants. Cluster galaxies have been classified into disks and spheroids 
on the basis of the shape of the light profile parametrized by the Sersic index n. The 
contamination between the two classes is expected to vary from ~ 8% at I = 22 to 
~ 17% at I = 22.5. 

Colour - magnitude distributions. Spheroids and disks show an evident segregation in 
the CM diagrams. The first family of galaxies define a sharp red sequence both in 
the optical and in the NIR, while disks have in general bluer colours and are located 
preferably below the CM relations. We find, however, some level of mixing between 
the properties of the two classes: some spheroids have blue colours, while some disks 
are red or redder than the red sequences. Interestingly, some of the blue spheroids are 
found at bright magnitudes (I < 20) and are not expected to be field contaminants. 
Spatial distributions. Disks and spheroids show a sharp spatial segregation, with disk 
galaxies found preferably in the low density regions. The spheroid population defines 
a central structure with a diameter of ~ 1.0 Mpc, and secondary clumps located at 
a distance >1.0 Mpc from the main density peak. By analyzing the density map of 
blue and red galaxies, we find that the two populations have a very different spatial 
distribution, the blue galaxies showing an amorphous structure. 

A standard Butcher - Oemler analysis has been performed in order to derive the frac¬ 
tion of blue cluster galaxies. The cluster EIS 0048 is characterized by a concentration 
parameter C ~ 0.45, that is intermediate between a very rich cluster and a less con¬ 
centrated, poor structure. The radius that encloses 30% of the total galaxy population 
is comprised between 1 and 1.5 arcmin (0.4 — 0.6 Mpc). The fraction of blue galaxies 
amounts to 0.11 ± 0.07, and increases up to ~ 23% at a distance of ~ 1.6 Mpc. We do 
not find, therefore, a strong BO effect in the cluster EIS 0048 at z = 0.64, in agreement 
with previous studies who did not find a high increase of the fraction of blue galaxies at 
z > 0.5. 
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